Bandwidth smearing in optical interferometry: 
Analytic model of the transition to the double fringe packet. 

R. Lachaume a ' 6 and J. -P. Berger c 

a Centro de Astroingenierfa, Departamento de Astronomfa y Astroffsica, Pontificia Universidad Catolica de Chile 

^ Max-Planck-Institut fur Radioastronomie 
c European Southern Observatory 

To appear in the proceedings of the 2012 SPIE Conference "Astronomical Telescopes and Instrumentation" 

ABSTRACT 

Bandwidth smearing is a chromatic aberration due to the finite frequency bandwidth. In long-baseline optical 
interferometry terms, it is when the angular extension of the source is greater than the coherence length of the 
interferogram. As a consequence, separated parts of the source will contribute to fringe packets that are not 
fully overlapping; it is a transition from the classical interferometric regime to a double or multiple fringe packet. 
While studied in radio interferometry, there has been little work on the matter in the optical, where observables 
are measured and derived in a different manner, and are more strongly impacted by the turbulent atmosphere. 
We provide here the formalism and a set of usable equations to model and correct for the impact of smearing 
on the fringe contrast and phase, with the case of multiple stellar systems in mind. The atmosphere is briefly 
modeled and discussed. 
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1. INTRODUCTION 

Long-baseline interferometry is based on the link between the complex visibility — the contrast and shift of the 
fringes — to the object's image, a Fourier transformEE Despite this straightforward relationship, the production 
of an image from optical interferometric data is difficult, even if a large number of measurements are recordedP 
One of the culprits is the turbulent atmosphere: It introduces a random variable optical delay, called piston, 
at each aperture, which in turn modify the phase of the complex visibility. At radio wavelengths, the problem 
has been overcome either by phase referencing on a close-by point-like sourceEE Also self-calibratiorSH can be 
used to retrieve the phases with the constraint they they cancel in a sum over telescope triplets^ — The sum 
of the phases over the triplet is called the closure phase. In the optical and the IR, the anisoplanatic patch 
in which the atmospheric phase can be considered uniform is too small to perform phase-referencing for most 
targets with the current sensitivity^ and the relatively small number of telescopes (typically 3-4) make the self- 
calibration difficult. For this reason, optical interferometrists still predominantly fit the visibility amplitude and 
closure phases with parametric models^ with imaging being relatively recent and based on blind reconstruction 
t echniques 11211111 

A problem arises when interferometers are used to observe wide fields, i.e. fields that are significantly larger 
than their point-spread function. Because of the finite bandwidth, the van Cittert-Zernicke theorem that links 
the visibility to the image is no longer exactly true and image reconstructions based on the fringes are smeared 
in a way that off-axis sources are elongated and less contrasted.^ As seen from an optical interferometrist's 
point of view, this bandwidth smearing occurs because (i) the finite spectral resolution R creates a fringe packet 
of finite width RX; and (ii) the different parts of the source are separated enough, by angular distance 9, for 
their respective fringe packets to overlap imperfectly, or even create a double fringe packet — the separation is 
B9 where B is the projected baseline. This effect becomes significant for BO ~ f?A/5P^ For a typical broad 
band near-IR filter at hectometric baselines, this corresponds to 5-10 milliarcseconds of separation when a wide 
band filter is used without spectral resolution. The use of shorter baselines can overcome the effect without a 
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significant precision loss on the separated parts of the source, yet this is not easily achievable nor desirable in 
practice, e.g. for hierarchical systems or in surveys for companions. 

While this bandwidth smearing has been well modeled in radio interferometry it is still relatively difficult to 
overcome a posteriori in image analysis.^ It is generally avoided using a higher spectral resolution. Alternatively, 
a posteriori software recombination with different phase centres using the amplitude and phase of the signal 
recorded at each antenna can solve the problem. At optical and IR wavelengths, higher resolution is generally 
barred by sensitivity issues and the nature of detection prevents to record the phase of the signal for later 
recombination, so the effect must be tackled. Despite the detailed description and modeling in the radio, the 
manner the visibility and related quantities are measured in optical and IR interferometry prevents us to use 
Bridle & Schwab's formalism^ as is. Since complex extended sources haven't been observed so far because of 
optical interferometers' limitations, we are concerned by the "simpler" multiple systems of near-point sources, 
that can be modeled using the visibility and closure phases. First steps were taken by Zhao et aiP^ who give an 
approximate formula for the smearing of the visibility amplitude of a binary and numerically model that on the 
closure phase. However, the derivation of generic formulas and a detailed analysis of the effect have not been 
performed yet. 

The purpose of this paper is to provide usable formulate for the smeared observables, square visibility am- 
plitude and closure phase, and assess the amount of bias that arises from this smearing in a practical case. We 
choose PIONIERP^ at the Very Large Telescope Interferometer. The instrument performs pairwise coaxial re- 
combination and delivers the interference pattern coded in temporal scans. We finally quickly assess the influence 
of the atmosphere using an analytic approximation of the turbulence. 

2. ANALYTIC DESCRIPTION OF THE FRINGE PACKET 

A point-like source o with flux density n*(£) at wave number a = (Tq + £ is observed by an interferometer 
in a spectral channel centred on a®. We consider the injection and coupling-corrected coherent flux, i.e. an 
interferogram whose continuum has been removed, the instrumental contrast of the fringes has been corrected 
for, and the imbalance of the flux coming from both arms has been removed. For coaxial recombination with 
temporal scans, it can be written as the sum over wavelengths and sources of point-source, monochromatic 
interferograms : 

M ab {5 ab {t)) = V / s a (C)?(0no(0e 2i,r(<T0+e)(x " i ' +5al,( * )+7r<,i,(t)) do-, (1) 

where 

• s a (£) is the complex spectral transmission of the electric field vector of incoming radiation in the interfer- 
ometer's arm, including terms that can be considered constant during a scan (< Is): quantum efficiency 
of the detector, absorption by the atmosphere, absorption along the arm, and instrumental optical path 
difference (OPD), i.e. instrumental phase. 

• Tr ab (t) is the atmospheric differential piston between arms a and 6, a term that varies during a scan (at the 
scale of ms) . 

• S ab (t) is the OPD modulation on baseline ab. 

• n o(0 i s the flux density of source o. 

• xf is OPD position of the centre of the fringe packets of object o. It is determined with the projected 
baseline p ab and object relative position in the sky using x ab = p ab ■ 6 . 

It can be rewritten as 

M ab {S ab {t)) = V s a ^n*(xf + S ab + ir ab (t)) e i< b +^™o(s ab +^"(t))^ 
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Figure 1. The telescope triplet and binary given in our test case. The telescope triplet a, b, c features baseline lengths B, 
0.6B, and 0.4B along the orientation of the binary, of flux ratio 0.6. 



• a"' = 2iraoXo b is the fringe packet shift in terms of phase. 

• / is the inverse Fourier transform of function /. 

3. PISTON-FREE OBSERVABLES 

In that case the coherent flux and its Fourier transform, the spectral density, read: 



M ab (5 ab ) = J2 s a s b n*(xf + 5 ab ) e ia ° b+2i ™ oSa 

O 

S ab {a)=Y,s a ^ b K(a-a )e 2 



J2iTtcrx n 



(3) 
(4) 



It is possible to measure s a s b on fringes obtained on a calibration lamp, provided that its spectrum is known, 
because piston is almost zero. We did it in the case of PIONIER, an example is shown in Fig. [2j Over the 
primary lobe of the fringe packet envelope the main trends seen on Fig. [2] are: 

• A constant instrumental phase over the primary lobe of the fringe packet (compare the flat parts of the 
curves in bottom right panel to the fringe envelope in the upper right one). 

• A clear group delay trend, with fringe packets centered at different positions in the different channels (top 
right panel). It is produced by the chromaticity of the instrumental differential phase (bottom left panel). 

In the following, we will approximate this by using a Gaussian profile of the spectral transmission, with a 
group-delay x°, b d and constant instrumental phase ip ab : 



s a s b {C) = exp 



s a s b (5) = exp 
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(5) 
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Figure 2. Example of transmission for PIONIER. The internal lamp was observed on 2011-08-09 in the H band with the 
3-channel spectral resolution, a sample baseline is reported here. The left columns display the spectral transmission and 
phase (s a s b ) for the three channels (from larger to lower wavelength: black & solid, red & dashed, and blue & dashed- 
dotted lines). The right columns shows the envelope and phase of the fringe packet (the Fourier transform of the latter). 
The phase slope of the spectral response (bottom left) translates into a shift of fringe packet centering, known as group 
delay (top right). 
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(b) Visibility amplitude as a function of binary separation 



200 




-200 I 1 1 1 1- 

5 10 15 20 



Reduced separation [mas • hm • /im~ 1 ] 

(c) Closure phase as a function of binary separation 

Figure 3. Example of a binary observed with an ideal interferometer of infinite spectral resolution (gray solid line) 
and with the three main spectral resolutions of PIONIER in the H band (7 channels, dark solid line; 3, gray dashed 
line; and whole band, dark dark line). The lines are easier to distinguish by level of smearing, the further from the 
ideal solid curve, the smaller the resolution. For the visibility and phase, the reduced binary separation is expressed 
in milliarcseconds-hectometres per micron for the sake of generality, i.e. in milliarcseconds for a hectometric projected 
baseline at a wavelength of one micron. For the closure phase the triplet of telescopes is taken with coordinates (0, B, 
0.4B) along the binary separation axis and the largest baseline B is used in the horizontal axis. In the whole band setting, 
the group delay is zero; in other settings we assumed a spectral channel with (-10, 0, 10) /im group delay. It corresponds 
to our PIONIER test case, shown in Fig. [2] 



Within this approximation, we take the example of a binary with flux ratio 0.6, and show its interferogram in 



Fig. 3(a) for the three spectral resolutions of PIONIER in the H band (1, 3, 7 channels, or 1Z ps 5, 15, 35) as well 
as the ideal interferometer are given in Fig. 3(a) See Fig. [I] for a sketch of the interferometer's configuration. 



3.1 Visibility 

The visibility amplitude can be determined using 

1 



|F|2 = iVW J R Mab ( S ) Mab ( § ) dS =N(vL Sat ^ Sab ^ d<7 ( ? ) 
where is a convenient normalisation factor. From Eq. Q, 

\V\ 2 = E ~ O (8) 

We now consider the case of the Gaussian transmission, assuming sources with a flat spectrum n*. The 
group-delay term of Eq. |5]) cancels out in the visibility amplitude formula so that: 

° °<P S_ - 

smearing 

where 1Z = ao/Aa is the spectral resolution. 

Within this approximation, the visibility of a binary with flux ratio 0.6 (see Fig. [I] baseline ab) for the three 
spectral resolutions of PIONIER in the H band (1, 3, 7 channels, or 1Z rj 5, 15, 35) as well as the ideal visibility 
are given in Fig. |3(b)| 

3.2 Closure phase 

The closure phase is the argument of the bispectrum, determined in direct space using 

abc / 7i fab ( cab /+\ \ /\ j be {xbci±\ \ njca 1 zca ( a\\ I kab /<\ the /i\ irca 



B abc = j Mf(5 a \t))M°^6 b ^t))M™{6 ca {t))\5 a \t)5 bc (t)S ca (t)\dt, (10) 

In the case of linear OPD variations we write S ab (t) ~ S ab t conveniently choosing the origin of time. It follows 
from Eqs. Jil&fio} and closure relation S ab + S bc + S ca = that 



B abc = J2 e K«f+<+0 f s a s b n*(xf + 6 ab t)s b ~^n*(x b p c + 5 bc t).^^ q {x c q a + 5 ca t) dt. (11) 

o,p,q ^ R 

It can be identified as a triple correlation and equated to the two-dimensional Fourier transform of a product: 

<xJ2<« (< b ~ &r>< - P a x c q a ) eK«S*+<+<*r). (12) 
with the triple spectral response 

(kabf: I kbet \ 
joe J' ( 13 ) 

In the case of a Gaussian throughput with linear group-delay (Eq. [5]) and a flat source spectrum: 

B abc x ^ - NoNpNge 1 ^ +<+<) e re^i°g2^^H^+^^ ^ (14) 

Smearing 

where ip abc = ip + ip bc + tp ca is the instrumental closure phase and af' = af — 1noQX ab A is the group- 
delay-corrected position of the fringe system for object o on base ab. It is worth noting that even with small 
separations (a° l « 0) the smearing term is non zero because of the group delay. Our example of a binary is 
display in Fig. |3(c) 



4. MODELLING THE ATMOSPHERE 



The extraction of the visibility amplitude and closure phase are meant to be robust against atmospheric tur- 
bulence, that is, they should not be biased. However, this has only been studied for compact sources where no 
smearing is present. We here quickly assess the combined influence of atmosphere and smearing. 

An analytic approach to the atmospheric turbulence can be taken in the slow turbulence regime, using the 
assumption that scanning is fast enough for the piston to vary linearly during one scan, i.e. ir ab = 7Tq 6 + TT^ b 5 ab , 
where 7rg h is the group-delay tracking error and ir^ b a rate of piston variation during scan. 7Tq 6 and 7r° b are 
random variables of zero mean. 

Using this approach, 

M ab (5 ab ) = J2 s^nfix^ + < b + (1 + ?rf )5 ab ) e ^M< b +^S b )+^o(i+< b )S a \ ( 15 ) 

o VI/ 

4.1 Visibility. 

The determination of the visibility is straightforward. The coherent norm stays identical, while the power 
spectrum in Fourier space involves a contraction of factor 1 + 7r" fc of Sf" b . Thus the only difference is a factor 
independent of separation (and smearing): 

|V 2 |jitter = TT ^|U 2 | id cal (17) 
1 + 7T° 

As long as irf b <C 1, the linear approximation ensures that it cancels out on a large number of scans. 

4.2 Closure phase in direct space. 

The piston terms cancel over a baseline triplet, so that the only changes concerns the triple correlation of the 
spectral response. The constant piston term can be seen as a fringe shift that adds to the original fringe shift of 
the object and the group delay 

x o ' x g.d. ^ x o + ^g.d. ' ~0 

and the linear variation of the piston can be seen as a scanning velocity change 

S ab ->■ S ab + nf 



The jittered bispectrum can thus be obtained from the piston-free formulas (Eqs. [T2 13 & 14) by performing 
these substitutions. 

It can be shown that a concave smearing function — as it is the case for the exponential in our approximation 
Eq. [14] — will tend to give a lesser value (i.e. more smearing) when averaged over zero mean random variables 
Trf b . We therefore expect that the atmosphere increases the smearing effects where they are already present. 
This is both confirmed in a work done by Zhao et alP^I with the Infrared Optical Telescope Array (but with a 
different data reduction method, thus bispectrum formula) and by a full numerical simulation of the atmosphere 
and interferometer that we are developing. 



5. CONCLUSION 

The smearing of the visibility and the closure phase derived from temporal scans can be modeled from a spectral 
calibration of the interferometer. We gave exact, generic analytic formulas and an application to a Gaussian 
band pass. This can be applied to "wide" multiple systems. When the components of the system are resolved 
analytic formulas can also be derived though they are slightly more complicated and do not bring insight into 
the phenomenon of smearing itself. 

The turbulent atmosphere has an impact on the amount of smearing of closure phase, though. While a 
reasonable analytic approximation can be done for each scan, there is no obvious method of taking atmosphere 
into account other than its average properties. Averaging the smeared observables over scans in the presence of 
the atmosphere needs some numerical modeling of the random atmosphere, which we are now undertaking. 



APPENDIX A. SOME QUESTIONS RAISED DURING THE CONFERENCE 



Is there an simple explanation of where the smearing comes from? It comes from the finite frequency 
band width. With infinite spectral resolution there is no smearing because fringes packets are infinitely long. To 
understand the phenomenon in optical interferometric terms, one can consider the extreme case of the double 
fringe packet: the fringes of one component are diluted by the incoherent flux of the other component, and the 
other way round, which leads to a visibility loss. 

What differs from the radio? We are measuring a square visibility amplitude instead of a complex visibility; 
the observables can be extracted from the interferogram in different, non equivalent manners; the atmosphere 
varies significantly during a scan. Formalisms look very similary, but there are significant differences in the 
formulae. 



Does it work for point sources only? It's possible to have compact resolved sources (i.e. each source, taken 
separately, is unsmeared) with little additional formalism. Also, sums can be substituted by integrals for more 
generality. 

Is it model dependent? Yes, you have to know what you source look like to model the smearing with this 
method. Also, you should know the spectrum of all parts of the source to accurately determine the fringe packet 
envelopes with accuracy; in the near IR, assuming a flat spectrum is a good enough approximation though (at 
least for young stellar objects). 

Why are there zeroes of the atmospheric bias in the smearing of the closure phase? I couldn't 
figure out the reason of this finding by Zhao et alP^ in the analytic formulae. 
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